Hamiltonian structure of classical iV-body systems 
of finite-size particles subject to EM interactions 

Claudio Cremaschini 
International School for Advanced Studies (SISSA) and INFN, Trieste, Italy and 
Consortium for Magnetofluid Dynamics, University of Trieste, Italy 

Massimo Tessarotto 
Department of Mathematics and Informatics, University of Trieste, Italy and 
Consortium for Magnetofluid Dynamics, University of Trieste, Italy 
(Dated: January 10, 2012) 

An open issue in classical relativistic mechanics is the consistent treatment of the dynamics of 
classical A-body systems of mutually-interacting particles. This refers, in particular, to charged 
particles subject to EM interactions, including both binary and self interactions (EM-interacting 
N-body systems). The correct solution to the question represents an overriding prerequisite for 
the consistency between classical and quantum mechanics. In this paper it is shown that such a 
description can be consistently obtained in the context of classical electrodynamics, for the case of 
a A-body system of classical finite-size charged particles. A variational formulation of the problem 
is presented, based on the A-body hybrid synchronous Hamilton variational principle. Covariant 
Lagrangian and Hamiltonian equations of motion for the dynamics of the interacting A-body system 
are derived, which are proved to be delay-type ODEs. Then, a representation in both standard 
Lagrangian and Hamiltonian forms is proved to hold, the latter expressed by means of classical 
Poisson Brackets. The theory developed retains both the covariance with respect to the Lorentz 
group and the exact Hamiltonian structure of the problem, which is shown to be intrinsically non- 
local. Different applications of the theory are investigated. The first one concerns the development 
of a suitable Hamiltonian approximation of the exact equations that retains finite delay-time effects 
characteristic of the binary and self EM interactions. Second, basic consequences concerning the 
validity of Dirac generator formalism are pointed out, with particular reference to the instant-form 
representation of Poincare generators. Finally, a discussion is presented both on the validity and 
possible extension of the Dirac generator formalism as well as the failure of the so-called Currie 
"no-interaction" theorem for the non-local Hamiltonian system considered here. 

I. INTRODUCTION 

In classical physics the formulation of the Hamiltonian mechanics of A^-body systems composed of interacting 
particles is still incomplete. This includes, in particular, the case of charged particles acted on by an externally- 
prescribed EM field as well as binary and self EM interactions. Indeed, based on general relativity (or special relativity, 
as appropriate in the case of a fiat Minkowski space-time) as well as quantum mechanics, common prerequisites for a 
dynamical theory for such systems should be: 

Prerequisite #1 : its covariance with respect to arbitrary local coordinate transformations. In the context of special 
relativity this requirement reduces to the condition of covariance with respect to the Lorentz group. 

Prerequisite #2: the inclusion of both retarded and local interactions. 

Prerequisite #3: the consistency with the Einstein causality principle. 

Prerequisite #4 '■ t ne validity of the Hamilton variational principle, yielding a set of equations of motion for all the 
A~ particles of the A^-body system. 

Prerequisite #5: the existence of a Hamiltonian structure. 

As clarified below, all of these statements should be regarded as intrinsic properties of classical A^-body systems 
which are characterized by non-local, i.e., retarded causal interactions, like those associated with EM fields [l|. In 
particular, requirements #4 and #5 involve the assumptions that the equations of motion of a generic A^-body system 
of this type should admit both Lagrangian and Hamiltonian variational formulations, obtained by means of a Hamilton 
variational principle, as well as a standard Hamiltonian form, i.e., a set {x, Hjy} with the following properties: 

a) x = (x«,i = 1, A^) is a super- abundant canonical state, with xW denoting an appropriate i-th particle canonical 
state; 

b) Hjy (to be referred to as system Hamiltonian) is a suitably regular function. In view of prerequisites #2 and 
#3, we expect to be prescribed in terms of a non-local phase-function of the form -Hat(x, [x]), x and [x] denoting 
respectively local and non-local dependences; 

c) for all particles i = 1,N belonging to the A^-body system, the variational equations of motion must admit 
the standard Hamiltonian form expressed in terms of the Poisson brackets with respect to the system Hamiltonian, 
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namely: 

^=x«,#J. (1) 

Here the notation is standard. Thus, (s(i), •••S(iv)) an d [t)i£\ = b7j£]( x ) are respectively the particles proper times and 
the local Poisson brackets (PBs). The latter arc defined in terms of the super-abundant canonical state x as 



<9x J — \dx 

with all components of x to be considered independent (i.e., x as unconstrained). Furthermore, J is the canonical 
Poisson matrix Q, while r)(x) and £(x) denote two arbitrary smooth phase- functions. 

Evidently, the above prerequisites should be regarded as overriding conditions for the transition from classical to 
quantum theory of the iV-body dynamics to be possible. However, despite notable efforts (see for example Dirac, 
1949 [3() the solution to the problem of fitting them together has remained still incomplete to date, at least in the 
case of systems of charged particles subject to EM interactions. 

From the point of view of classical physics the reason is related to the nature of EM interactions occurring in 
A^-body systems. These can be carried respectively both by external sources (unary interactions, due to prescribed 
external EM fields) and by the particles themselves of the system (internal interactions). The latter include both 
binary EM interactions acting between any two arbitrary charged particles and the self EM interaction, usually known 
as the EM radiation- reaction (RR; Dirac Q, Pauli @, Feynman Q). As earlier pointed out (see Refs.Q and 0], 
hereon referred to as Paper I and Paper II respectively) a rigorous treatment of the EM self-interaction consistent 
with prerequisites #l-#5 can only be achieved for extended classical particles, i.e., particles characterized by mass 
and charge distributions with finite support. In particular, a convenient mathematical model is obtained by assuming 
that these classical particles are non-rotating and their mass and charge distributions are quasi-rigid in their rest 
frames (see Paper I for a detailed discussion on this point). In Papers I and II the dynamics of single extended 
particles (1-body problem) in the presence of their EM self-fields was systematically investigated in the context of 
classical electrodynamics, by means of a variational approach based on a Hamilton variational principle. As a result, 
the Hamiltonian description for isolated particles subject to the combined action of the external and the self EM 
interaction has been established. 

However, fundamental issues still remain unanswered regarding the analogous formulation of a consistent dynamical 
theory holding for classical A^-body systems of finite-size charges subject to only EM interactions (EM-interacting 
N-body systems). In fact, it is well known that traditional formulations of the relativistic dynamics of classical charged 
particles are unsatisfactory, at least because of the following main reasons. 

The first one is related to the approximations usually adopted in classical electrodynamics for the treatments of 
RR phenomena. In most of previous literature charged particles are regarded as point-like and the so-called short 
delay-time ordering 

e = — « 1 (3) 

is assumed to hold, with t and t' denoting respectively the "present" and "retarded" coordinate times, both defined 
with respect to a suitable Laboratory frame. This motivates the introduction of asymptotic approximations, both for 
the EM self 4-potential and the corresponding self-force, which are based on power-series expansions in terms of the 
dimensionless parameter e and are performed in a neighborhood of the present coordinate time t (see related discussion 
in Paper II). Nevertheless, previous approaches of this type have lead in the past to intrinsically non-variational and 
therefore non-Hamiltonian equations of motion These are exemplified by the well-known LAD and LL RR 

equations, due respectively to Lorentz, Abraham and Dirac (Lorentz, 1985 10]; Abraham, 1905 [ll[ and Dirac Q) 
and Landau and Lifschitz Both features make these treatments incompatible with the physical prerequisites 

indicated above. 

The second motivation arises in reference to the so-called "no- interaction" theorem proposed by Currie [l3l [l4j . 
According to its claim, isolated A^-body systems, formed by at least two point particles, which are subject to mutual 
EM binary interactions and in which the canonical coordinates are identified with the space parts of the particle 
position 4-vectors, cannot define a Hamiltonian system with manifestly covariant canonica l eq uations of motion. The 
correctness of such a statement has been long questioned (see for example Fronsdal, 1971 [15| and Komar, 1978-1979 
[TB43)- in particular, the interesting question has been posed whether the "no-interaction" theorem can actually be 
eluded in physically realizable classical systems. Interestingly, Currie approach is based on the well-known generator 
formalism formulated originally by Dirac (DGF; Dirac, 1949 Q). Therefore, related preliminary questions concern 
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the conditions of validity of DGF itself and, in particular, whether such an approach actually applies at all to EM- 
interacting N-body systems. 

Goals of the paper 

Put all the previous motivations in perspective, in this paper a systematic solution to these issues is presented, 
for the case of EM-interacting TV-body systems of classical finite-size charged particles. The theory is developed in 
the framework of classical electrodynamics and special relativity (i.e., assuming a flat Minkowski space-time) and is 
shown to satisfy all the prerequisites indicated above (#l-#5). 

Starting point is the determination of both Lagrangian and Hamiltonian equations of motion for classical charged 
particles belonging to an EM-interacting A^-body system. By construction (see Section 2) the latter can be considered 
as a system of smooth hard sphere, namely in which hard collisions occurring between the particles, when their 
boundaries dQ^) (see below) come into contact, conserve each particle angular momentum. In particular, for simplicity, 
in the following all extended particles will be considered as acted upon only by EM interactions, thus ignoring the 
effect of hard collisions on the A^-body dynamics. As in previous Papers I and II, where the 1-body problem was 
investigated, the derivation is based on the variational formulation of the problem. This requires, in particular, 
the determination of the appropriate variational functionals required for the description of both binary and self EM 
interactions. The resulting action functional is found to be expressed as a line integral in terms of a suitable non- 
local variational Lagrangian, the non-locality being associated both to the finite-extension of the charge distributions 
and to delay-time effects arising in binary EM interactions. Based on the Hamilton variational principle expressed 
in superabundant variables, the resulting variational equations of motion (for the A^-body system) are then proved 
to be necessarily delay-type ODEs. As a consequence, based on the definition of suitable effective Lagrangian and 
Hamiltonian functions, a manifestly-covariant representation of these equations in both standard Lagrangian and 
Hamiltonian forms is reached. The main goal of the paper is then to show that these features allow a Hamiltonian 
structure {x, Hjy} to be properly defined in terms of a suitable superabundant canonical state x and a non-local 
system Hamiltonian function Hjy. The result follows by noting that the Hamiltonian equations in standard form 
admit also a representation in terms of local PBs, defined with respect to the super-abundant canonical state x. 

A remarkable development concerns the construction of an approximation of the Hamilton equations in standard 
form. This holds in validity of both the short delay-time and large-distance orderings, namely under the same 
asymptotic conditions usually invoked in the literature for the asymptotic treatment of the RR problem. Based on the 
analogous approach developed in Paper II, it is shown that a suitable N-body Hamiltonian approximation [of the exact 
problem] can actually be reached, which preserves its Hamiltonian structure. In particular it is proved that, unlike in 
the LAD and LL equations, the asymptotic approximation obtained in this way keeps the variational character of the 
exact theory, retaining the standard Lagrangian and Hamiltonian forms of the A^-body dynamical equations as well 
as the delay-time contributions arising from the various EM interactions. 

Further interesting conclusions are drawn concerning the validity of DGF in the present context. This refers, in 
particular, to the so-called instant-form representation of Poincare generators for infinitesimal transformations of the 
inhomogeneous Lorentz group. It is pointed out that DGF in its original formulation only applies to local Hamiltonian 
systems and therefore is inapplicable to (the treatment of) the EM-interacting A^-body systems considered here. 
For definiteness, the correct set of Poincare generators, corresponding to the exact non-local Hamiltonian structure 
{x, Hn} determined here, together with their instant-form representation, are also provided. This permits to develop 
a modified formulation of DGF, denoted as non-local generator formalism, which overcomes the previous limitations 
and is applicable also to the treatment of non-local Hamiltonians in terms of essential (i.e., constrained) canonical 
variables. 

Finally, on the same ground, the Currie "no-interaction" theorem is proved to be violated in any case by the 
Hamiltonian structure, i.e., both by its exact realization {x, Hpj} and its asymptotic approximation. Counter-examples 
which overcome the limitations stated by the "no-interaction" theorem are explicitly provided. In particular, the 
purpose of this discussion is to prove that indeed a standard Hamiltonian formulation for the Af-body system of 
EM-mutually-interacting charged particles can be consistently obtained. The main cause of the failure of the Currie 
theorem is identified in the conditions of validity of DGF on which the proof of the theorem itself is based. 



Scheme of the presentation 

The paper is organized as follows. In Section 2 the expressions for the EM current density and self 4-potential for 
extended charged particles are recalled. In Section 3 the A~-body action integral is constructed, based on the results 
of Papers I and II and focusing in particular on the derivation of the action integral for the EM binary interactions. 
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Section 4 deals with the dynamical equations of motion in Lagrangian form for the interacting particles, derived 
from the TV-body hybrid synchronous Hamilton variational principle (THM.l). A standard Lagrangian form of the 
same equations is obtained, by introducing an effective non-local Lagrangian function (Corollary to THM.l). In 
Section 5 the corresponding non-local Hamiltonian theory is developed in terms of a variational Hamiltonian function 
(THM.2). A standard Hamiltonian form of the TV-body dynamical equations is constructed and the associated non- 
local Hamiltonian structure is proved to exist (Corollary to THM.2). General implications of the non-local TV-body 
theory are presented in Section 6, while in Section 7 the TV-body Hamiltonian asymptotic approximation of the exact 
solution is developed (THM.3). In Section 8 the validity of DGF is investigated, in connection with the non-local 
Hamiltonian structure determined here and the corresponding instant-form representation of Poincare generators 
implied by it. The extension of DGF to non-local system is presented in Section 9, by means of the introduction of 
the notion of non-local Poisson brackets. Then, an overview of the basic features of the "no-interaction" theorem is 
given in Section 10. Explicit counter-examples to the theorem are determined, which are provided by the existence 
of standard Hamiltonian forms for locally-isolated 1-body systems and globally-isolated TV-body systems (THM.4). 
The failure of the "no-interaction" theorem and the non-applicability of its statements to the non-local Hamiltonian 
TV-body system considered here are discussed in Section 11. Finally, Section 12 provides a closing summary of the 
main results, while the mathematical details concerning the derivation of the TV-body action integral for the binary 
interaction are reported in the Appendix. 



II. TV-BODY EM CURRENT DENSITY AND SELF 4-POTENTIAL 



For the sake of clarity, let us first briefly recall the key points of the formulation presented in Papers I and II 
dealing with the definitions of extended charged particle and of the corresponding 4-current and self 4-potential. We 
shall assume that each finite-size particle is characterized by a positive constant rest mass ttiq and a non- vanishing 
constant charge gW, for i = 1,TV, both distributed on the same support dQu\ (particle boundary). More precisely, 
for the i-th particle, the mass and charge distributions can be defined as follows. Assuming that initially in a time 
interval [— oo, t ] the i-th particle is at rest with respect to an inertial frame (particle rest-frame 1Z where the external 
forces acting on the particle vanish identically), we shall assume that: 

1) In this frame there exists a point, hereafter referred to as center of symmetry (COS), whose position 4-vector 
r (i)COS = ( c ^ r o) spans the Minkowski space-time A4 4 C R 4 with metric tensor i)^ =diag(+l — 1 — 1 — 1). With 

respect to the COS the support dfi^) is a stationary spherical surface of radius > of equation (r — r c ) 2 = er^. 

2) The i-th particle is quasi-rigid, i.e., its mass and charge distributions are stationary and spherically-symmetric 
on dQ(i) (concerning the physical requirements assuring the condition of rigidity we refer to the related discussion in 
Paper I). 

3) Mass and charge densities do not possess pure spatial rotations. Therefore, introducing for each particle the 
Euler angles a(s(j)) = V'}( S( . ) ) which define its spatial orientation (see definitions in Ref. 20]), the condition of 
vanishing spatial rotation is obtained imposing that a(sm) = const, holds identically. As a consequence, only the 
translational motion of charged particles need to be taken into account. 

In addition, as stated before, hard collisions occurring between the particles are considered ignorable. As a conse- 
quence, for all particles the equations of motion (flj are assumed to hold identically, namely for all s^ S I = M. 

For each extended particle the covariant expressions for the corresponding charge and mass current densities readily 
follow (see Paper I). In particular, these can be expressed in integral form respectively as: 

,(0, r+°° 



(r) 



dsu Mll (s)5(\x(i ) \ - o- (i ))5(s - si (l )), 



(4) 



3 { ±s(r) = tV / dsu^(s)S(\x (i) \-a (i) )S(s-s 1{i) ), 

47rfJ (i) J-oo 



where by definition Siu\ is the root of the algebraic equation 



Ml* 



(»i(i)) 



= 0. 



Here the notations are analogous to those given in Paper I, so that in particular !i"' i (s(i)) 
of the COS for the i-th particle, while 



x (i) -r r 



( s w) 



(5) 



(6) 



is the 4- velocity 



(7) 
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Finally, following the equivalent derivations given in Papers I and II, the non-divergent EM self 4-potential 
Aj^ el ^^\r) for the single (namely, i-th) extended particle can be readily obtained as well. It is sufficient to re- 
port here the solution for A^ el ^^ (r) which is valid in the external domain with respect to the spherical shell of the 
same particle. For a generic displacement 4- vector Xwv £ M 4 of the form 

X ^ = r » - r®» (s W ) , (8) 

which is subject to the constraint 

xW"uW( flw )=0, (9) 

this sub-domain is defined by the inequality 

X^XjpK-ofq. (10) 

In such a set, A^ el ^ % ' (r) is expressed in integral form by the equation 

r 2 ^ 

A( sel W)(r) =2g« / dr'5{R^ a R$), (11) 



where RW a is the bi-vector 

flW« = r o_ r W» s (12) 

with r^ % ' ,a = r^' a (s'^) being the i-th particle COS 4- vector evaluated at the retarded proper time s',^, obtained as 
the causal root of the equation R^ a Ra = 0. As remarked in Paper II, Eq.lfTTj) formally coincides with the analogous 

solution holding for a point particle. However, in difference with the point-particle case, A^ e (r) can be defined 
everywhere in M 4 in such a way to be non-divergent (see Paper I) . 



III. THE NON-LOCAL iV-BODY ACTION INTEGRAL 



In this section we formulate the TV-body Hamilton action functional suitable for the variational treatment of a 
system of N finite-size charged particles subject to external, binary and self EM interactions. In such a case, in 
analogy to Paper II, the action integral can be conveniently expressed in hybrid super-abundant variables as follows: 

S N (r,u,[r})= [^(^«) + 4 eXt)W M+4 Sei/)W (^M) + 4 bm)W (^W)] (13) 

i—l.N 

(non-local action integral). As in Papers I and II, r and it represent local dependences with respect to the 4- vector 
position and the 4-velocity, while [r] stands for non-local dependences with respect to the 4-vector position. In 
particular, the latter are included only via the functionals produced by the EM-coupling with the self and binary 
EM fields for the i-th particle, namely S^ eL ^^ and 5^ OT )W. Instead, S^} and Sq^^' identify for each particle the 
functionals produced by the inertial mass and by the EM-coupling with the external EM field. We stress that the 
functionals S^(r,u), S^ xt ^ l \r) and S^ el ^^' (r, [r]) are formally analogous to the case of a 1-body problem treated 
in Papers I and II and can be represented as line- integrals (see below) . We now proceed evaluating explicitly the new 
contribution S^! m ^ l \r, [r]). 



A. 5^'™"'' (r, [r]): EM coupling with the binary-interaction field 

The action integral g^ m )M ( rj [ r ]) containing the coupling between the EM field generated by particle j, for j = 1,N, 
and the electric 4-current of particle i is of critical importance. Its evaluation is similar to that of the action integral 
of the self-interaction outlined in Paper I. For the sake of clarity, in this subsection we present the relevant results, 
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while the details of the mathematical derivation are reported in the Appendix. According to the standard approach 
[HI, (r, [r]) is defined as the 4-scalar 

4 Wn)(<) (*',M)= £ 4" m)fo) (r,[r]), (14) 

3=1,N 

where (r, [r]) is defined as 

4 Wn)W) (^. M) = f dn- 2 A^ l W{r)W (r) , (15) 

with j4( sei ^^^ M (r) being the EM 4-potential generated by particle i at 4-position r, whose expression is given by 

Eq. (fTT|) . In addition, jj;P (r) is the 4-current carried by particle j evaluated at the same 4-position and given by 
Eq.Q, while dtt is the invariant 4- volume element. In particular, in an inertial frame Si with Minkowski metric 
tensor r/^, this can be represented as dtt = cdtdxdydz, where (x,y,z) are orthogonal Cartesian coordinates. As 
shown in the Appendix, an explicit evaluation of the action integral (|15l) yields the following representation: 

S$ in) M(r,[r}) = 2 J^lj\ r y ( S(i) )^ 2 ^>( S(j) )5(^>^) -af,)), (16) 

where and syj are respectively the proper times of particles i and j, while RS l i' a denotes 

fl(tf)a = r (i)a ( S{ .^ _ r ««( S(i) ). (17) 

It is worth pointing out the following basic properties of the functional S^ n,t '^'. First, it is a non-local functional 
in the sense that it contains a coupling between the "past" and the "future" of the particles of the iV-body system. 
In fact it can be equivalently represented as 

4«"><«>(,,h) . «wl r «, b ^m r^i^*^-^,. da) 

c J — oo J — oo U) 

Furthermore, the iV-body system functional (|14l) is symmetric, namely it fulfills the property 

J2 S^\r A ,[r B ])= £ S^\r B ,[r A ]), (19) 

i,j = l,N i,j=l,N 

where ta and tb are two arbitrary curves of the iV-body system. 

B. The non-local iV-body variational Lagrangian 

In this section we provide a line-integral representation of the Hamilton functional Sn in the form 

/+oo /-+oo 
ds (i) 4\r,u,[r])= ]T / T (i) (r,«,[r]), (20) 
%=i,jn -°° i=l,N J -°° 

where T^(r, [r] ,u) and L^(r, [r] ,u) are respectively the i-th particle non-local contributions to the fundamental 
Lagrangian differential form and to the corresponding non-local variational Lagrangian. Invoking Eq. (|18|) and recalling 

also the results of Paper II, (r, [r] , u) can be written as 

L?{r,u, [r]) = L$(r,u) + L^ t)W (r) + L^ /)W (r, [r]) + 4 &m)W (r, [r]), (21) 

where L^(r,u), L^ xt '^'(r) and L^ eL ^^ l \r, [r]), x^ m ^^ (r, [r]) denote respectively the local and non-local terms. In 
particular, the first one is the contribution carried by the iner tial term, while L^ m , L%* m and L^ m)(l) identify 
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respectively the external, self and binary EM-field-coupling Lagrangians. These are defined as follows: 

~dr^ 1 



L^(r,u) ee m«cu» 



ds(i) 2 



J c ( r ) = — -r — A u (r w (s(i)),o- (i) ), 



(r, [r]) 



Hi) 

q(i) dA 1 ^ —(self)(i) 

-A,, , 



c (is 



(0 



(bin) (i) 
C 



(r, M ) - £ 4 MW) (Mr]) = ^^ £ ^r )W) '- 



C 067,") 



(22) 
(23) 
(24) 
(25) 



Here, A*"* ^\ A^" 1 ^^ an( j ^( 6m )W) d eno t e the surface-averages performed on the i-th particle boundary 9fi(i) 

(see Paper I) respectively of the external, self and binary EM 4-potentials. In particular, j[^ el ^^ anc ] ^( 6m )( y ) are 
defined as 



.4 



(self)(i) 



2g« / drf'5{R^R^ - a 2 {{) ), 



ds 



dr»(s (j) ) 



ds 



In addition, is the bi-vector 



(26) 
(27) 

(28) 



with S(j) and sC^ denoting respectively "present" and "retarded" proper times of the i-th particle. 



IV. NON-LOCAL TV-BODY VARIATIONAL PRINCIPLE AND STANDARD LAGRANGIAN FORM 

Let us now proceed constructing the explicit form of the TV-body relativistic equations of motion for each extended 
charged particle in the presence of EM interactions (i.e., including external, binary and self EM interactions). This 
is achieved by adopting for the TV-body problem a synchronous variational principle [2l|, [22j which, in analogy with 
the approach developed in Papers I and II, can be expressed in terms of the super- abundant hybrid (i.e., generally 
non-Lagrangian) variables 



(29) 



and for a suitable functional class of variations {/}. The latter is identified with the set of real functions of class 
C fc (R), with k > 2, and fixed endpoints which are prescribed for each particle i = 1,TV at suitable proper times su\i 
and S(i) 2 , with s (i)1 < s (i)2 , i.e., 



{/} 



/ (i) (^)i) = /f; 

i = l,iV; j = 1,2 and k> 2 



(30) 



It follows that by construction the variational derivatives of the Hamilton functional Sn (see Eg . (1201 ) are performed 
in terms of synchronous variations, i.e., by keeping constant the i-th particle proper time suy The result is expressed 
by the following theorem. 

THM.l - TV-body hybrid synchronous Hamilton variational principle 

Given validity of the prerequisites #l-#5 for the N-body system, let us assume that: 



1. The Hamilton action SN(r,u, [r]) is defined by Eg. \20\) . 
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2. The real functions f^(su\) in the functional class {/} [see Ea. (30p ] are identified with the super- abundant 
variables i29\) which are subject to synchronous variations Sf^(su\) = /^( s (i)) — /i ( s (i))- The l a tt er belong 
to the functional class of synchronous variations {Sf^}, with 



(31) 



(32) 



for fc = l,2,V/W( S(l) ),A W ( S(l) )G{/}. 

3. The extremal curves /^(sfi)) € {/} for Sn, which are solutions of the E-L equations 

6S N (r,u, [r]) = 
*/W(*(0) 

exist for arbitrary variations Sf^(s^) (hybrid synchronous Hamilton variational principle). 

4. If the curves r^ M (s(i)) ; /or £ = 1, JV are a/Z extremal, each line element ds^ satisfies the constraint 

dsl^n^dr^is^dr^is^). 

5. The 4-vector field A^f (r) is suitably smooth in the whole Minkowski space-time AT 1 . 

6. The E-L equations for the extremal curves r^^(s^) are determined consistently with the Einstein causality 
principle. 

7. All the synchronous variations Sfj^ (sm) (k=l,2 and i = 1,N) are considered as being independent. 

It then follows that the E-L equations for and r^' M following from the synchronous hybrid Hamilton variational 
principle i32\) give respectively 



(33) 



' ' V = m^cdr^ - m^cu^^dsu) = 0, 

r III ° ° \ l l ' 

du 



(0 

SS 



N 



7« 



oY«M( S(i) ) 



mWcd«W( Sw ) + ^ifj>*)(O dr (0»( Si ) = 0, 



where is the total Faraday tensor acting on particle i and given by 



p(tot)(i) ,-p( sel f){i) ~p{bin)(i) 



(34) 
(35) 

(36) 



where all quantities are intended as surface-averages on the i-th particle shell-surface dQu\. Eqs. \3J$ and i35\) are 
hereon referred to as N-body equations of motion. In particular: 

1) F ' (r^'j = d^A v — d v A^ is the antisymmetric Faraday tensor of the external EM field evaluated on 

the extremal curve = rWP • 

F ( r ^\ [ r ^]) = Fju/ 1 ^ (r^ (s(i)) ( s (i))) * s ^ e non-local antisymmetric Faraday tensor produced 
by the EM self- field of the i-th particle and acting on the same particle. This is given by 



pM/)(0 



namely 



F 



(self)(i) 



2q 



^(s'^'-u^js'^ 



(0 „,«, 



5(0 



(0 



(37) 



(38) 



s ( i)= s (i)- s <i)"t 



where the delay-time S(j\ ret is the positive (causal) root of the 1-particle delay-time equation 



(39) 
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^ui/ * s ^ e non-local antisymmetric Faraday tensor produced on particle i by the action of all the remaining 



particles, i.e. F^ n) %) = T,j=l,N F ^ n)ilj) (r-W , [r^] , ct (i) , cr (i) ) , w/iere 



,cr (i); cr (i) 



.ffere </«e notation is as follows. Huu is defined as 



s U)= s (ij)\' T U)) 



7 (i) 



i?fe>n^(s (i) ) 



(is 



0') 



(40) 



(41) 



while the delay-time S(j) — si^ ( <7 (i)) an d s^j) — s^fl «?~e respectively the positive (causal) roots of the 2-particle 

delay-time equations 



R {i)« R {i) 



.2 

(<) 
2 

(J) 



0. 
0. 



(42) 
(43) 



Therefore, sj^ and s|^j depend respectively on and <ryj . 

Proof - The proof is analogous to the corresponding 1-body problem detailed in THM.l of Paper 1. Indeed, since 
the Dirac-deltas S^R^^R^ — a?*) and 5{R^^ a Ra^ — cr?--,) are independent of particle 4-velocities, the variations 



with respect to u£ deliver necessarily the E-L equation (|3"4"1) . To prove also Eq. ([33)) . we notice that the synchronous 
variations of the functionals S^{r,u), S^ xt ^ l \r) and S^ el ^^\r, [r]) necessarily coincide with those of the 1-body 
problem (see Papers I and II). Therefore, it is sufficient to inspect the variational derivative of the non-local binary- 
interaction functional 5^"" (r, [r]). Its variation with respect to Sr^^^s^) takes the form 



5S. 



(bin) (i) 
C 



J2 {[SA + 5B] m + [SA + SB} m } 



(44) 



j=l,N 



where 



5B {IJ) ee S^tfe^drC^^ 



(45) 



and the second term [5 A + 5B\^^ follows by exchanging the particle indices. Then, using the chain rule and integrating 
by parts, after elementary algebra Eqs. (|40|) and (|41l) follow (for details see Appendix A in Paper I). In agreement 
with the Einstein causality principle the positive roots of the delay-time equations (|4"2"j) and (|4"3"|) must be selected. 
Q.E.D. 

A few comments are here in order regarding the implications of THM.l. 

1. Coordinate-time parametrization of the N-body equations of motion 

It is important to stress that for each i-th particle, its equations of motion [in particular the E-L Eqs.(f34"|) and 
(|35l) ] can be parametrized in terms of the single coordinate time t rather than the corresponding particle proper 
time S(j). This is obtained introducing the representations in terms of the single coordinate (i.e., Laboratory) 
time t e I e K, namely letting for all i = 1, N 



-Mm 



it) 



ds 



(ct,rW), 
cdt 

7(i)' 



with 7^) and f3^ denoting the usual relativistic factors 



7(0 = 



-1/2 



)/c. 



(46) 



(47) 
(48) 
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This implies also that the 4-velocity can be represented as u^" 1 
equations f[34|) and (|35|) become respectively 

m^cdr®* - m^cv (i)ti dt = 0, 



h(i) vil) ' 1 with 



-rn 



7 W 



(c, vW). Hence, 

(49) 
(50) 



2. Delay-time effects 

Delay-time effects which appear both in the EM RR and binary interactions are due to the extended size of the 
charged particles. In particular, the delay-time characterizing the self-interaction acting on particle i depends 
only on the radius of the charge distribution of the same particle. Instead, the delay-time appearing in the 
binary interaction experienced by particle i depends either on the radius am of particle i or on the radii o~u-\ of 
all the remaining particles. In the case of TV-body system of like particles, such that o~u-\ = au\ — a, the two 
terms on the r.h.s. of Eq. (|44[) coincide yielding a single delay-time contribution in Eq. (|40[) . Explicit evaluation of 
delay times involves the construction of the positive (causal) roots of the equations (|42l) and (|43| . Based on the 
coordinate-time parametrization (|46|) . these can be solved explicitly for the coordinate delay-time t^ ret = t',* —t. 
The causal roots are in the two cases respectively 

*(*)«*(*) = V [r«(t)-rW(t - * Wret (t))] 2 + ^ > 0, (51) 



t (ij)ret (t) = - c yj [rW(t)- r W)(i - t (lj)ret {t))} 2 + af }) > 0. (52) 

Notice that the same roots can also be equivalently represented in terms of the corresponding particle proper 
times (s/j) for i = 1,N). For this purpose it is sufficient to introduce for the i-th particle proper time the 
parametrization = su\(t) which is determined in terms of the coordinate time t by means of the equations 
([55)1 . It follows, in particular, that the proper delay-times corresponding to (|5~T1) and (|5"2"|) become respectively 

S(i)ret{s(i)) = s(t(i)ret(t)} an d s {ij)ret{s(j)) = s(*(y)ret (t)). 



Along the lines of the approach given in Paper II, it is immediate to show that the hybrid- variable variational 
principle given in THM.l can be given an equivalent Lagrangian formulation. An elementary consequence is provided 
by the following proposition. 

Corollary to THM.l - Standard Lagrangian form of the iV-body equations of motion 

Given validity of THM.l, let us introduce the non-local real function 

Le fi ,N= Y, £ i//(W r D. (53) 

i=l,N 

where L e ff t x is denoted as N-body effective Lagrangian, while nfl^{r,u, [r]) is defined as 

L%(r, u, [r]) ee lW(r,«) + L^\r) + 2L%*™{r, [r]) + L$™{r, [r]). (54) 
Here L$, I^ xt '^' and [Jg^f'w coincide with the variational Lagrangians defined above (see Eas. \22\l - \2J$ ) , while 

T (bin)(i) ■ . , 

L e ff is given by 

L^irM) - E 2qit)qU) """I M r^^ 1 ^, (55) 
n c ds (i) J-oo ds u) 

with K^i' being the sum of Dirac-deltas 

KM ee 5{R^ a R^ ] - <J 2 U) ) + 5(R<M a RM - o- 2 (i) ). (56) 
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Then, the E-L equations {3J$ and H35\) coincide with the E-L equations in standard form (see Papers I and II) 
determined in terms of the N-body effective Lagrangian L e tf t jf. In particular, the E-L equations in standard form 
for the i-th particle become 

9Leff ' N -0, (57) 



F^\r)L eff , N = Q, (58) 



where 



F®( r ) = — - - (59) 

<**(<) 

denotes the E-L differential operator. 

Proof - The proof is based on THM.2 of Paper I and by noting that the E-L differential operator acts only on local 
quantities. Hence, the equivalence of Eqs. ([57|) -([58 | with the corresponding E-L equations (f34|) and (|35|) follows by 



elementary algebra and in view of the identity F^ (r)L e ff n — (r)L *Jf. 
Q.E.D. 

To conclude this Section a final remark is necessary regarding the functional setting of the iV-body equations of 
motion given above. 

We first notice that the E-L equations (|34p and ([55]) , and the equivalent Lagrange equations in standard form (|58l) , 
imply for alH = 1,N the second-order delay- type ODEs 

mfc = S_ F (totm, dr (^ {Si y m 

ds (i) c 

Let us introduce the Lagrangian state w = ((rW,u^) ,i = 1,JV) , with = — , spanning the iV-body phase- 
space Tn = Yl -^1(0) where Tui) = M^f) x U$ and U$ = K 4 indicate respectively the corresponding 1-body phase 

i=l.N 

and velocity spaces, the latter endowed with a metric tensor n^ u . To define the initial conditions, let us make use 
of the coordinate-time parametrization (|46|) . denoting w(t) = ((f^^ (t), u'^ (t)) ,i = 1,N) and r^'(t) = r^\s^(t)), 
u^'(t) = «W(sjj)(f)). Then, a well-posed problem for Eqs. (|60p is obtained prescribing the initial history set {w} t(j C 
L n ■ For an arbitrary coordinate initial time to G I = R, this is defined as the ensemble of initial states 

{w} fo ={((f«(t), ««(*)) eC^- 1 \l),i = l,N),yte [*o-CT(*o).*o].*>2}. (61) 

Here, for a given initial (coordinate) time to G I, i™* x (io) denotes the maximum (for all particles) of the delay-times 
t(i)ret(*o) and t(y) re t(*o)) namely 

i£f (to) = max{< Wret (to),t w(ret (to), Vi, j = l,iV} . (62) 

Solutions of Eqs. (|5D|) fulfilling the initial conditions defined by the history set {w} tQ are sought in the functional class 
of smooth 4-vector solutions of the form r^' 1 = r^^s^)), with su\ = su\(t) and t € Iq = (to, 00 ) > which belong to 
the functional class 

{r( s )} = { (r«",i = l,iv)| r«" = r«"(fl W ) £ C«(7), with S(i) = 8(i) (t) e (7^(1), fc > 2,Vi e /„} • (63) 

In the following we shall assume that m £/ie setting defined by Eq. A63\) with the history set i61\), the ODEs A60\) admit 
a unique global solution of class C^ k ~ x \lo) , with k > 2. 



V. A?-BODY NON-LOCAL HAMILTONIAN THEORY 



Based on THM.l and its Corollary, an equivalent non-local Hamiltonian formulation can be given for the hybrid and 
Lagrangian- variable approaches stated in THM.l and Corollary. The strategy is similar to that developed in Paper II. 
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Thus, first we proceed constructing the intermediate set of hybrid variables y = (r,p) = (r^ipft = l,iv) and the 
related non-local variational Hamiltonian = i/{; l '(r,p, [r]), identified, as usual, with the Legendre transformation 



(i) 

of the corresponding non-local variational Lagrangian L\ . Hence, for all i = 1,N: 



" als 



1 ! 



while pffl is the i-th particle conjugate momentum defined in terms of as 



P W = ^1 



ds (i) 



From THM.l it follows that 



(i) 



pW - m^cuf + 



(64) 



(65) 



(66) 



where A^ ot ^^ is given by 



-r(6m)(y) 



according to the definitions given in Eas.(|2"6" )) -([2"7 |) . As a consequence, H\' becomes simply 



j=l,N 
(i) 



H®(r,p, [r]) = 



2ml i) c 



„(»') _ L.i(«)(i) 



(')m _ £ ^(tot)(i)/i 



This permits us to represent the Hamilton action functional in terms of the hybrid state y, yielding 

SH N (r,p, [r])= S hP> 

i=l,N 



(67) 



(68) 



(69) 



where 



S H u) = 



S(i)2 



ds 



(i) 



M ds {l) 



(70) 



represents the ?-th particle contribution. Of course, as an alternative, analogous dynamical variables can be defined 
also in terms of the effective Lagrangian L^lt [see Eq. (f54|) ]. This yields the notion of effective Hamiltonian H^h 

and of the corresponding state x = (r,P) = (r^ 1 ^ ,P^\i = l,N^j , which will be shown below to identify a (super- 



abundant) canonical state (see Corollary to THM.2). Thus, H^h - to be considered a non-local function of the form 

l eff ~ 2J eff\ 



— H^j '^(r, P, [r]) - is prescribed in terms of the Legendre transformation with respect to L^Jf, namely letting: 



H (i) _ p(i) 



e " ~ M ds 



J eff> 



(71) 



while denotes the effective canonical momentum 



p{i) = 

dsta 



(72) 



From the Corollary to THM.l it follows immediately that 



7 (i) 



(73) 
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where A[ to l\[ 1 ^ is given by 



Atot)(i) _ -r{ext)(i) -j{self){i) 
(e//)(U A* A* 



— (6in)(ij) 



(74) 



and 



-j-(bm)(ij) 



+ oo 



ds 



dr»{s {j) ) 



ds 



(./) 



(75) 



with ifW) being given by Eq.(|56p. Finally, H^jf becomes 



H 



(i) 



1 



Therefore, by direct comparison with Eq. (|68|) it follows identically that 



p(i)fJ. _ 1 4( tot )( i )A t 



(76) 



H 



Then the following theorem, casting the Hamilton variational principle of THM.l in terms of the state x, holds. 



(77) 



THM.2 - iV-body non-local Hamiltonian variational principle 

Given validity of THM.l with Corollary and the definitions nW - t7(% ) as well as (7Jp - |73p , let us assume that the 
curves /^(s(,)) = yW = (rW^jty s )( S( .j) belong to the functional class {/} of C 2 -functions subject to the boundary 
conditions 



for k = 1, 2, S(£)xj s (i)2 S/CR and wii/i S(£)i < S( 2 ;)2- T/iera i/ie following proposition holds: 
The modified Hamilton variational principle 



SS 







(78) 



(79) 



subject to independent synchronous variations <5/W(s(j)) = (8r^^{s^),Sp^\s^) [ j performed in the functional class 
indicated above, yields, for all i = 1,N, the E-L equations 



5S H (i) 



(80) 



5S h™ 



= 0. 



(81) 



These equations coincide identically with the N -body variational equations of motion \3Jj.\) and \35\) . Hence, the set 
{y,-ff/v} = \y^ l \H\\i — lj-^l defines a non-local Hamiltonian system. 

Proof - The proof is analogous to that given in THM.2 of Paper II. In particular, it can be reached, after elementary 
algebra, by invoking the symmetry properties of the variational functional Sh n > namely 



Sh n (ta,P, M) = S HN (r B ,P, [ta]), 



(82) 



where again and rs are two iV-body arbitrary curves of the functional class {y}. It follows that the variational 
derivative in the E-L equation Eq. (|5Tj) becomes 



SS 



SS H^ 



5S„ y 
8 [r^H 



(83) 
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where the summation is performed only on the non-local contributions. As a consequence, the E-L equations (|80[) 
and (ED) yield 



SS H^ 

Spy 



ds 



(i) 



pf 



c A » 



= 0, 



ss 



dp^ 



ds 



o /1 (tot)(i) 



Q r {i)v 



= 0. 



(84) 
(85) 



Taking into account the definitions given by Eq. (f66)) the equivalence with Eqs.([34| and (|35|) is immediate. 
Q.E.D. 

Let us now pose the problem of the construction of the corresponding TV-body Hamiltonian equations in standard 
form, as suggested by the results of Paper II. The non-local Hamiltonian system {y, H^} is said to admit a standard 

Hamiltonian form |x, H^ffi ■■■^^ffi ^ the -/V-body equations of motion can be cast, for all i = 1,N, in the form 



dr (l ^ dH 



0) 

eff 



ds 



dP, 



(<) ' 



(86) 



dP, 



(0 



dH 



(0 

eff 



ds 



(i) 



(87) 



(i) - 

in terms of a suitably-defined effective particle Hamiltonian H e Jf, to be identified with Eq. (|76|) . In particular, a 
./V-body system with state x = (x",i = l,iV} is said to be endowed with a Hamiltonian structure {x, iJjv, e //} if, for 
all particles belonging to the iV-body system, the equations of motion for the i-th canonical particle state can be 
represented in the PBs notation ^ in terms of a single Hamiltonian function i?/v,e//, i-e., for all i = 1,N 



ds 



(i) 



N,eff 



with HN te ff denoting a still to be determined, appropriate system effective Hamiltonian. Extending the treatment 
holding for the 1-body problem (see Paper II), here we intend to prove that the Hamiltonian structure {x, iJ/v>//} 
holds also in the case of EM-interacting TV-body systems. The following proposition holds. 

Corollary to THM.2 - Standard Hamiltonian form and Hamiltonian structure of the iV-body equa- 
tions of motion 

Given validity of THM.2 and the definitions given by Eqs. J76p, it follows that: 

TC2i) The non-local Hamiltonian system {x, Hpj} admits a standard Hamiltonian form defined in terms of the 
set |x, ...H^jfJ \, with H^ff the i-th particle effective Hamiltonian [given by Eq. (|76[) ]; furthermore x = 
{x( ,; \i = l,iV} is the super- abundant canonical state, while 



(rW°,r«) 



x« 



(89) 
(90) 
(91) 

are respectively the i-th particle canonical state, ^-position and effective canonical momentum [defined by Eq. (|73l) ] . 
As a consequence, Eqs. A80\) and &81\) can be cast in the standard Hamiltonian form A86}) and J&7] ). 

TC22) The equations \86\) and admit also the equivalent representation \88\) and hence the set {x, iJ/v,e//} 
defines a Hamiltonian structure, with iijv.e// being the effective N-body Hamiltonian function 



H 



N.eff 



H 



(i) 
eff 



(92) 
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TC23) Introducing the system Hamiltonian 



H 



N 



i=l,N 



(0 



(93) 



defined in terms of the variational i-th particle variational Hamiltonian H^' [see Eq.(|&S|)]. it follows identically that 



Hn = H 



N,eff- 



(94) 



Proof - TC2i) The proof follows from straightforward algebra. The first equation manifestly reproduces Eg. (1501) . 
because of the definition of fl^y given above. Similarly, in the second equation the partial derivative of H^y, recovers 

TC22) To prove the existence of the Hamiltonian 



the correct form of the total EM force expressed in terms of ^['//j^ 



OH 

structure, it is sufficient to notice that — j£y- = [r^^, iJ/v.e//] and 



0H { 



3P> 



drWf 



N,eff 



TC2 3 ) By construction 

[see Eqs.(|68| and (|76p ] for alH = 1,N the effective and variational Hamiltonians coincide [see Eq. (f77)) ]. This implies 
the validity of Eq. ([94| too, namely Hm identifies also the system Hamiltonian. It follows that for all particles i — 1, N 
the canonical equations of motion (|88[) recover the standard Hamiltonian form expressed in terms of the PBs with 
respect to the the system Hamiltonian, i.e., Eqs.([T]), so that {x,i/jv} identifies the Hamiltonian structure of the 
EM-interacting TV-body system. 
Q.E.D. 



VI. GENERAL IMPLICATIONS OF THE NON-LOCAL W-BODY THEORY 

Let us now comment on the general implications of the previous theorems. 

Remark #1 - Difference form of the Hamilton equations of motion. The canonical equations ([I} imply the following 
difference equations, i.e., the infinitesimal canonical transformation generated by H^: 



dx w = ds 



(0 



N 



(95) 



Remark #2 - Coordinate-time representation of the Hamilton equations of motion. Introducing the coordinate-time 
parametrization ([4^)) . Eqs. (P5|) become 



cdt r 

7(i 



N 



(96) 



with 7^) denoting again the relativistic factor (|4"7|) . while [•, •] are the local PBs evaluated with respect to the super- 
abundant canonical state x. These yield explicitly 



dt 1 

7(i) m£ ] 



p(0 



dt 

7(i) 



caf=dtvf, 



(97) 



,u (i) /) / 4( to *)( i ) / (0 \ ji ( «4(*°*)( i ) a) «/i( to *)(0 

7(i) c dA 1 ^ c dr^ 1 * 



dpi. 

7(») TOo C ar 



-4 l 



Remark #3 - Well-posedness of the N-body equations of motion. All the equations of motion indicated above [see 
THMs.l and 2 and their Corollaries] are equivalent to each other and are manifestly Lorentz-covariant (see also 
related discussion in Paper II). Then, a well-posed problem for the Hamiltonian equations in standard form (|88[) can 
be obtained in analogy to the problem defined by Eqs. (160l) .(|61 [l . This is achieved, first, by prescribing the appropriate 
initial history set {x} to C IV. For an arbitrary coordinate initial time to G I = K, this is defined as the ensemble of 
initial states 



{x} to ee {S(t) ee ((?%),£«(*)) e C { - k - 1 \l),i= l,N) ,Vt € [to -Cr(*o)>*o],fc> 2} 



(99) 



16 



(canonical history set), where t^ x (to) is the maximum delay-time at to [see Eq. ([62|) ]. Furthermore, in analogy 
with Ea.(|55]). solutions of Eas. ([55|) fulfilling the initial conditions defined by the history set {x} to are sought in the 
functional class (|30]) by identifying 



(100) 



In the following we shall assume that in the setting defined by Eq. \10O\) with the canonical history set h99f) , the ODEs 
\88\) admit a unique global solution of class C^ _1 ^(7o) with k > 2. 

Remark #4 ~ Extremant and extremal curves. In all cases indicated above the solutions of TV-body E-L equations of 
motion (extremal curves) and of the Lagrangian and Hamiltonian equations in standard form given by the Corollaries 
to THMs.l and 2 (extremant curves), satisfy identically, for all i — 1,N, the kinematic constraints 



ttJT'uW" = i 



(velocity constraints) and 



(101) 



(102) 



(line-element constraints). The first constraint implies that the time-components of the 4- velocity depend on the 
corresponding space components, while the second constraint requires that the particles' proper times are uniquely 
related to the corresponding coordinate times. In particular, we shall denote as extremant canonical curves 



x(s (1) ,...s w ) = {x (1) (s (1) ),...x (Ar) (s (Ar) )j 



with = xW(s/j)) for all i = 1, N, arbitrary particular solutions of the canonical equations (|123[) . 
Remark #5 - Unconstrained varied functions. By assumption, both the varied functions /W = 



and X W = (r^,P^) ( 



(103) 



(««) 



entering respectively THMs.l and 2 as well as the Corollary of 



•'"' = "■ ' " ; - >( S(l) ) auu = ^•"■>- f 7i ';(a (<) ) 
THM.2 arc unconstrained, namely they are solely subject to the requirement that end points and boundary values 
are kept fixed (and therefore do not fulfill the previous kinematic constraints). This implies, in particular, that all of 
the (8) components of /W, yM and x^ must be considered independent. On the other hand, both the extremal and 
extremant curves satisfy all of the required kinematic constraints, so that only (6) of them are actually independent 
for each particle (see also discussion in Paper I). 

Remark #6 - Non-local Hamiltonian structure and unconstrained canonical state. Thanks to proposition TC2 3 of 
the Corollary to THM.2 the Hamiltonian structure {x, Hjy^ff} coincides with {x,Hn} , Hn denoting the non-local 
system Hamiltonian defined by Eq. ([93f . We remark, however, that in the PBs given by Eq.fl| the partial derivatives 
must be evaluated with respect to the unconstrained states x^ and not indicated above. This means that Hjq 
must be considered a function of x. It is immediate to prove that the same Hamiltonian structure {x, Hn} holds 
provided the super-abundant canonical state x = (xS^,i — 1,N) is considered unconstrained. In fact, as shown 
by the Corollary to THM.2, in such a case the canonical equations in standard form (|86[) and ([87]) admit a PB- 
representation of the form (|88|) . For this purpose let us make use of the coordinate-time parametrization x = x(t), 
denoting x(t) — x Q + dx, with x Q = x(t a ) and dx = (dx^ 1 ', ...dxW) . Furthermore let us require that the initial 
history set {x} to is prescribed. Then, a necessary and sufficient condition for the equations of motion to admit the 
standard Hamiltonian form (JXJ) is that the fundamental local PBs for the state x = x(t), defined with respect to the 
same state x G = x(t a ), namely 



p(») pti) 
r ( i )M j pU) 



(x„) 

(Xo) 
(Xo) 



= o, 
= o, 



(104) 



are identically satisfied for all i,j = 1,N and \x,v — 0,3. This is realized only when the super-abundant variables 







ljA^) are considered independent. 



Remark #7 - The canonical flow is not a dynamical system. A final issue concerns the properties of the flow 
generated by the canonical problem (|88|> and (|M)) - (|100I) (canonical flow). In the iV-body phase-space this is an 
ensemble C < -' c ~ 1 '— homeomorphism (with k > 2) of the type 



(105) 
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which maps an arbitrary history set {x} to C IV onto a state x(£) crossed at a later coordinate time t (i.e., at t > to). 
This map does not generally define a dynamical system. In fact, unless there is a subset on non- vanishing measure of 
r^r in which {x} t(j reduces to the initial instant set 

{x} to = {x(t ) = x } , (106) 

the flow (|105[) is not a bijection in F^. To prove the statement it is sufficient to notice that - in the case of a non-local 
Hamiltonian structure {x, HN.eff} - if the history set is left unspecified and only the initial state x(£ ) is prescribed, 
the image of the initial state is obviously generally non-unique [and hence it may not coincide with x(i)]. In fact, 
while the same initial state x(i Q ) may be produced by different history sets, for example {x} to and {x'} to , the same 
history sets will generally give rise to different images x(£) and x'(i). 

We emphasize that for TV-body systems subject to EM interactions the instant set (|106[) can be realized only for 
special initial conditions, i.e., for example, if for all t < to all the particles of the system are in inertial motion with 
respect to an inertial Lorentz frame. Since, unlike the external EM field, binary and self EM interactions cannot be 
"turned off" , it follows that the set of initial conditions (|106[) has necessarily null measure in IV . 



VII. THE TV-BODY HAMILTONIAN ASYMPTOTIC APPROXIMATION 



In this section we want to develop asymptotic approximations for the equations of motion of EM-interacting TV-body 
systems. This involves different asymptotic conditions to be imposed on both the self and binary interactions. In 
particular: 

1) For the RR self-interaction of each particle i with itself this is provided by the short delay-time ordering, namely 
the requirement that the dimensionless parameters e^) = _fi!L_fM_ ) f or £ = 1, TV are all infinitesimal of the same order 
e, i.e. e ~ e^) <C 1, S(j) — s',^ denoting the i-th proper-time difference between observation (s) and emission (s') of 
self-radiation. 



2) For the binary EM interactions the Minkowski distance 



R(ij) 



between two arbitrary particles of the system is 



much larger than their radii, in the sense that for all i,j = 1, TV, with i ^ j, the large-distance ordering < 
< e holds. 



The fundamental issue arises whether an approximation can be found for the TV-body problem which: 

1) is consistent with the orderings 1) and 2); 

2) recovers the variational, Lagrangian and Hamiltonian character of the exact theory (see THMs.l and 2); 

3) preserves both the standard Lagrangian and Hamiltonian forms of the equations of motion; 

4) retains finite delay-time effects characteristics of both the RR and binary EM interactions, consistent with the 
prerequisites #l-#5 indicated above. 

In this regard, a fundamental result is the discovery pointed out in Paper II of an asymptotic Hamiltonian approx- 
imation of this type for single extended particles subject to the EM self-interaction. This refers to the retarded-time 
Taylor expansion of the Faraday tensor contribution carried by the RR self- force. More precisely, in the case of a 
single particle, this is obtained by Taylor-expanding the RR self-force 

0?-^ir ,m (r» W ,^(^))*^M (m 

for i = 1 (see Eq. (|38|l l in the neighborhood of the retarded proper-time s'^y Here we claim that an analogous conclusion 
can be drawn also for the corresponding TV-body problem, by introducing the same expansion to all charged particles 
and invoking the large-distance ordering for the binary interaction. For this purpose, we shall assume that the external 
force acting on each charged particle is slowly varying in the sense that, denoting r' = r' 4 ^ f s (j)J an( l r = r ^ >1 ( s (i))i 



( 



F { ;f(r') -f;;;i,i| , - om. nm 



c 



(r')-F% Xt) 


(r) r 


- 0(e), 


)-F { ;f{r) 




- 0(e), 


-F^fir)) 


,hk 


- 0(e). 
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Then, the following proposition holds. 

THM.3 - iV-body asymptotic Hamiltonian approximation. 

Given validity of THM.2 and the short delay-time and large- distance asymptotic orderings as well as the smooth- 
ness assumptions fil08}) - HTTW for the external EM field, neglecting corrections of order e™, with n > 1 (first-order 
approximation), the following results hold: 

T3\) The vector fields JiflTj ) describing the RR self- force are approximated in a neighborhood of as 




(111) 



to be referred to as retarded-time Hamiltonian approximation for the self- force, in which the first term on the r.h.s. 
identifies a retarded mass-correction term, rn^^j = J,^ - denoting the leading-order EM mass. Finally, g$ are the 
^-vectors 

«W<o))~^ 

T,?2^ The tensor fields F^ v ^ J for all i, j — 1, N , withi^j, appearing in the binary EM interaction ( see Eg. (JUty ) 
are approximated by the leading-order (point-particle) terms: 

,<r (i) = (VT (i) =0). (113) 

T3z) The corresponding asymptotic N-body equations of motion obtained replacing G$ and F^ n * 5 with the 
asymptotic approximations illl]) and \113]) are variational, Lagrangian and admit a standard Lagrangian form. De- 
noting with rip' = ro ^ s '(;)) ^ e extremal i-th particle world-line at the retarded proper time s'^, the i-th particle 
asymptotic variational Lagrangian functions become: 




-p(bin){ij) ^ —(bin){ij) 



Lts ym (r,u, [r]) = L«(r,«) + 4~*>W(r) + ) + L<£M{r, [r]). (114) 

Here L^J and L^ xt remain unchanged (see Eqs. (22]) and (23))). while the non-local terms "^-m ( r ' M) an< ^ 
^(>vm( r i M) are respectively 

^a'^VVf) = ff »(r«')r^, (115) 
T (bin)(i), rn q® dr®* ^ -r(«n)(«) / n N 

^C.J^M) = — L A M K")= )' ( 116 ) 

W j=l,JV 



where, from Eq. (27\ ) one obtains 



°(j) 



Similarly, the effective particle Lagrangians are, for i = l,N: 

^ha Sym (r, u, [r]) ee L«(r, u) + 4-*>«(r) + ^^VV^) + 2LgJ>» (r, [r]). (118) 

T3±) The N-body equations obtained imposing the asymptotic approximations given by Eos, kill}) and \113\) are 
also Hamiltonian. The asymptotic variational and effective Hamiltonian functions are given respectively by 

H W - D (0*^_iW (119) 
= pW_ tW fi20l 
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with L^ asym and defined by Eqs. {114-\j and H118\) , while now 



orCO 
(i) Uasym 



PW = . (122) 

Proof - T3i) The proof is analogous to that given in THM.5 of Paper II. T3 2 ) To prove the validity of Eq. (fTT3|) . 

let us recall the definition of ^ given by Eq. (|4"0)) . Then, for each particle, imposing the large-distance ordering 

and neglecting corrections of order e™, with n > 1, the leading-order contribution is given by Eq. (|113p , which depends 
on a single delay-time determined by the positive root of the equation R^ a Ra^ = 0. T33) The proof follows 
by first noting that the function 1^'^' contributes to the i-th particle E-L equations only in terms of the local 
dependence in terms of r^\ Second, in the large-distance ordering, the asymptotic approximation for the iV-body 
Lagrangian carrying the binary interactions yields a symmetric functional, as the exact one. Therefore, the N- 
body asymptotic equations are necessarily variational and Lagrangian. Straightforward algebra shows that the E-L 
equations determined with respect to the asymptotic variational Lagrangian (|114|) coincide with the asymptotic 
approximations proved by propositions T3i) and T3 2 ). In a similar way it is immediate to prove the validity of 
Eq. (|118| ). which shows that the same asymptotic equations admit a standard Lagrangian form. T34) Finally, the 
equivalent iV-body variational and standard Hamiltonian formulations follow by performing Legendre transformations 
on the corresponding asymptotic variational and effective Lagrangian functions. It follows that the asymptotic iV-body 

(i) 

equations of motion can also be represented in the standard Hamiltonian form in terms of f asym • 
Q.E.D. 



It is worth pointing out the unique features of THM.3. These are related, in particular, to the asymptotic expansion 
performed on the RR self-force alone. In most of the previous literature, the short delay-time expansion is performed 
with respect to the particle present proper-time. This leads unavoidably to local asymptotic equations (analogous to 
the LAD and LL equations) which are intrinsically non- variational and therefore non-Lagrangian and non-Hamiltonian. 
In contrast, the short delay-time expansion adopted here (as in Paper II) approximates the non-local RR vector field 
in a manner that meets the goals indicated at the beginning of the section. The remarkable consequence is that the 
asymptotic TV-body equations of motion retain the representation in standard Hamiltonian form characteristic of the 
corresponding exact equations. Finally, we stress that in all cases both for the exact and asymptotic formulations, 
the variational and effective Lagrangian and Hamiltonian functions are always non-local functions of the particle 
states. The non-locality is intrinsic and arises even in the 1-body systems, being due to the functional form of the 
EM 4-potential generated by each extended particle. 



VIII. ON THE VALIDITY OF DIRAC GENERATOR FORMALISM 



A seminal approach in rclativistic dynamics is the Dirac generator formalism developed originally by Dirac (Dirac, 
1949 Q) to describe the dynamics of interacting iV-body systems in the Minkowski space-time. Dirac's primary goal 
is actually to determine the underlying dynamical system, exclusively based on DGF. In his words ll In setting up such 
a new dynamical system one is faced at the outset by the two requirements of special relativity and of Hamiltonian 
equations of motion". It thus "... becomes a matter of great importance to set up (in this way) new dynamical systems 
and see if they will better describe the atomic world" (quoted from Ref.Q). 

DGF is couched on the Lie algebra of Poincare generators for classical iV-body systems. The basic hypothesis behind 
Dirac approach is that these systems must have a Hamiltonian structure {z, Kn} of some sort, with z = {z^ 1 -* , ...z^ } 
and Kn being a suitable canonical state and a Hamiltonian function of the system. In particular, the canonical states 
z^ of all particles i — 1, N must satisfy, by assumption, covariant Hamilton equations of motion of the form 



ds(i) 



A' 



(123) 



with s/j) denoting the z-th particle proper time. However, it must be stressed that certain aspects of Dirac theory 
remain "a priori" undetermined. This concerns the functional settings both of the canonical state z and of the 
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Hamiltonian function Km- In particular, the system state z = {z' 1 ', ...zW} remains in principle unspecified, so that 
it might be identified either with a set of super- abundant or essential canonical variables. Thus, for example, in the 
two cases the i-ihe particle state z^ l > might be prescribed respectively either as: a) the ensemble of two 4-vectors 
zW = M')' i ,7rW' i ) , with r"' 1 = W'' ,r''') an d 7r^^ = (7rW°,7rW) being respectively the particle 4-position and its 
conjugate 4-momentum; b) the ensemble of the corresponding two 3- vectors obtained taking only the space parts of 
the same 4-vectors and 7r^^, namely in terms of z'W = (rW , tj-W) . Thus, depending on the possible prescription, 
the very definitions of the PBs entering Eqs. (|123[) and DGF change. 

Furthermore, it remains "a priori" unspecified whether is actually intended as a local or a non-local function of 
the canonical state z. In Ref. 3], however, certain restrictions on the nature of the set {z, K^} are actually implied. 
These will be discussed below, leaving aside for the moment further discussions on this important issue. 

Provided the Hamiltonian structure {z,K^} exists, any set of smooth dynamical variables 77, £ and £ depending 
locally on the canonical state z necessarily fulfills the following laws 

[77, £] = - [£, 77] , 

fat + C] = ht} + h(], 
[€,vC} = l€,vK + v[Z,C], 

[[^],C] + [b?,CU] + [[C,£U]=o. (124) 

DGF relies on the Lie transformation formalism and is based on the representation of the Lorentz transformation 
group in terms of the corresponding generator algebra. This is defined as the set of phase-functions (Poincare algebra 
generators) {F} given by 



F = -v>"a lx + -M^b llv , (125) 



with a M , 6 M „ being suitable real constant infinitesimals and p*", — —M v ^ appropriate local phase- functions 

obeying the PBs (Lorentz conditions) 



[Pv,Pv] = 0, 



M^Ma/j = -r] tia M V + ri va M^ -r]^M au + r} u pM alit . (126) 

Hence, F as given by Eq. (|125l) generate respectively infinitesimal 4-translations (for b^ v = and ^ 0) and 4- 
rotations (for 6 M „ ^ and a M = 0, corresponding either to Lorentz-boosts or spatial rotations) via infinitesimal 
canonical transformations of the type 

z^z' + (5 z, (127) 

with S z ~ 0(S), 5 > denoting a suitable infinitesimal. 5 a z is determined identifying it with 

S z = [z, F] , (128) 

to be referred to as the local variation of z. Hence, according to DGF an arbitrary dynamical variable £ depending 
locally and smoothly on the canonical state z transforms in terms of the law 

£ e = £ + <«, (129) 
S ^z) = [^,F}, (130) 

where S !;(z) = [£(z') — £(z)] [1 + 0(8)] . It is immediate to determine an admissible representation for the generators 
{F} = j^ 1 ", . Let us first consider a relativistic 1-body system represented by a single particle in the absence 

of external forces. For definiteness, let us assume that Vs^) 6/e1 the particle 4-velocity is constant. Then p^ and 
M^&ie manifestly given by: 

Pv = rf\ 

= ^^W-g^TrW, (131) 

where respectively the 4-vectors qjp and are to be identified with the 1-body coordinates and momenta. In 
the case of the corresponding TV-body problem for relativistic interacting particles, in Dirac paper three different 
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realizations of were originally proposed, which are referred to as the instant, point and front forms. All 

of them follow by imposing the velocity kinematic constraints (|10ip . In particular, the instant form is realized by 

prescribing the reference frame in such a way to set 7q = 0, for all i = 1, N, namely describing each particle position 

only in terms of the space components rW = of its position 4-vector, for I = 1,3. In detail, recalling Eq. (l90l) and 
introducing the notation 

^ = (4\^% (132) 

according to Dirac the instant form (for TV-body systems of interacting particles) is obtained by imposing the velocity 
kinematic constraints (|101| on the free-particle canonical momenta 7rj^ ee M = m^cu^ = f7rjy ee0 , 7rj>eeJ P- e -' m * nc 
absence of an external EM field] such that 



4L,o = ^/m^ + 4 ) r l (133) 

and then introducing a suitable interaction ^-potential = (Vo,V) taking into account all the particle interac- 
tions. Letting l,m = 1,3, this yields the iV-body Dirac constrained instant-form generators (pQ,pi, Mi m , Nio^j 
represented in terms of the constrained states z'W = (rW,7i"W) (for i — 1,N): 



Po 



Pi 



E P8°= E +*$l + Vo, (134) 

i=l,JV i=l,N 



E *r. ( 135 ) 

i=l,N 

M lm = E [^-^M ]. ( 136 ) 



= 1,JV 



Mo = E W + «W+Vi, (137) 

i=l,N 

with Vi = Vq ^2a = i r / an< i denoting the time-component of a suitable interaction potential 4-vector = (Vo, V). 
Here both po and iV;o are still expressed in terms of the free-particle canonical momentum 7Tyy ee , while A/jjo differs 

from M;o because of the imposed kinematic constraint. Therefore, if interactions occur, their contribution show up 
only in po and Niq. In Ref. [|[ the interaction-dependent Poincare generators were called "Hamiltonians" . 

Nevertheless, for the validity of the transformation laws (|129p as well as of the Lorentz conditions (|126[) . Eqs. (|134[ )- 
(|137p are actually to be cast in terms of the 4-momenta of the interacting system 7fW (rather than the free particle 
momenta 7r2 ee )- This means that in general 7i"W should be considered as suitably-prescribed functions of 7rV? ee and 
of the interaction 4-potential V^. For definiteness, let us consider the case of an isolated TV-body system subject only 
to binary interactions occurring between point particles of the same system. In such a case the interaction potential 
4-vector is necessarily separable [2{|, i.e., such that 

v ll = (y ,v) = E v ¥ } > ( 138 ) 

i=l.N 

(i) . (i) 

with V/j, denoting the i-th particle interaction potential ^-vector. Then, assuming that are only position- 
dependent, in view of Eqs. (|134l) and (|137[) , for each particle the canonical 4-momentum of interacting particles ir^ 

(i) (i) 

must depend linearly on and 7r^ ree , namely it takes the form 



tJP=*52Um-^ { °. (139) 
which implies in turn that necessarily ir^ — J mh^c 2 4- (7rW — VM) 2 + V® . As a consequence, Eqs. (|134|) - (|137|) are 
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actually replaced with 

po 



E p [ o } = E v m ° )2c2 + - vW ) 2 + ( 14 °) 

i=l,N i=l,N 

% = E *f. ( 141 ) 



=1,JV 



M lm = 2 [^-rWTrW]. ( 142 ) 



i=l,JV 



- E r * (I) v m ° l)2c2 + - vW ) 2 + v ^ ( 143 ) 

i=l,JV 

where (| 141 [) and (|142[) retain their free-particle form. In order that the Poincare generators po and pi commute [in 
accordance with Eqs. (ll26[) . with PBs now defined in terms of the constrained state z'], then it follows necessarily that 

the 4- vectors V^, for all i = 1,N, must be local functions of the 4-positions of the particles of the TV-body system, 
namely 

V« = V«(rf 1J ,...rW). (144) 

The explicit proof of this statement is given below [see the subsequent subsection #3 and in particular the inequality 
(|155p ]. Hence, by construction, in the Dirac approach the "Hamiltonians" po and Niq are necessarily local functions 
too. It follows that DGF applies only to local Hamiltonian systems. 

It is worth noting that the same conclusion follows directly also from Dirac's claim (see the quote from his paper 
given above) that the iV-body system should generate a dynamical system. In fact, in the customary language of 
analytical mechanics the latter is intended as a parameter-dependent map of the phase-space r^v onto itself. This 
means that, when the canonical state z is parametrized in terms of the coordinate time t, there should exist a 
homcomorphism in r^r of the form: 

2(t )=*o++a(t), (145) 

with t, t a S J = K. Therefore, if the previous statement by Dirac is taken for granted, in view of the discussion 
reported above (see Remark #7 in Section 7), it implies again that Km must only be a local function of the system 
canonical state z. 

Let us now analyze, for comparison, the implications of the theory developed in the present work for EM-interacting 
TV-body systems. 



1. Non-local Hamiltonian structure 

The first issue is related to the Hamiltonian structure {z,Kn} which characterizes these systems. According to 
the Corollary of THM.2 this should be identified with {x, Hn} . Thus, the super-abundant canonical state z should 
coincide with x = {xW,i = l,iV} spanning the 87V-dimensional phase-space Tn = Hi=i t N^i (with c M 8 ), 
while Kn is identified with the non-local Hamiltonian Hn = Hn{t, P, [r]). In particular, x^ 1 ^ = (r^^jPjfij is 

the corresponding i-th particle canonical state, with and denoting respectively its position and canonical 
momentum 4-vectors. In the present case and in contrast to DGF, it follows that: 

1. Property #1: the system Hamiltonian Kn = Hn must be necessarily non-local. 

2. Property #2: the super-abundant state x = (x^ l ',i — 1,N\ is canonical, namely it satisfies the canonical 
equations JT|) in terms of the local PBs defined with respect to the same state. This occurs if x is considered 

unconstrained, i.e., when the i-th particle state xW is identified with xW = P^^j . Hence, the Hamiltonian 

structure {x, Hjy} holds in the unconstrained 8iV-dimensional phase-space Tn = ^ii=i.N^i \ with c M 8 . 
As a consequence, also the PBs (including the fundamental PBs (| 104|) and the Lorentz conditions (|126[1 ) are 
defined with respect to the unconstrained state x. 

3. Property #3: only the extremal or extremant canonical curves x(sm, ..S(jy)) [see Eq. (|103[) ] and not the varied 
functions satisfy identically the kinematic constraints f|101[) and (|102[) . 
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4. Property #4 '■ the 4-potential must be necessarily a non-local function. In particular, for binary EM interac- 
tions it must be a separable function, i.e., of the form (|138|) : 



V,{rM) = E VjpirM), (146) 



i=l,N 

W«(r, W ) = ^4efft (147) 



c 



with Affijd being defined by Eq.<[73]l. 



2. Conditions of validity of Dirac instant-form generators 



A further issue is related to the representation of the Poincare generators and, in particular, to the instant form 
representation given by Dirac and usually adopted in the literature. The latter is based on Eqs. (|134|) - (|137|) . rather 
than on Eqs. (|140l) - (|143|) . in which tt^ 2 replaces 7r2g e under the square root on the r.h.s. of Eqs. (ll34[) and (|137|) . On 
the other hand, based on the non-local Hamiltonian structure {x, Hjy} expressed in terms of the unconstrained super- 
abundant canonical state x [see Eqs. ([5§|) . (|90D and (|9ip], an admissible realization for can be determined 
which holds for an arbitrary N > 1. In fact, it is immediate to verify that the phase- functions 



E p , 



(i) 



i=l,N 

%» = E [rfW-Wlf 

i=l.N 



(148) 



satisfy identically the PBs (|126l) expressed in terms of the same state x. In view of Property #2 this requires that 
the canonical generators defined by Eq. (|148p must be considered independent. Hence, no constraints (on 

them) can possibly arise by imposing the validity of the PBs (|126[) . 

Another possibility, however, lies in the adoption of a constrained formulation. This is obtained imposing the 
kinematic constraints f|101[) and identifying the canonical state with the constrained vector x' = (x'W,i = l,iV) with 
X 'M — (rW,PW). Recalling again Eqs. ([90|) and (f9Tj) . here rW and denote respectively the space parts of the 
corresponding i-th particle 4- vectors. 

To carry out a detailed comparison with Dirac, let us consider in particular the instant-form representation of 

as given by Eq. (|148l) . In such a case the generators are represented by the set, defined for l,m = 1,3: 



Pa 
Pi 

Mlm 

Nio 



E^ 

i=l, N 



E *P> 

= 1,N 
=1,JV 

E H^o ( 



' m 'ml 



(149) 
(150) 
(151) 
(152) 



i=l,JV 



where poj Pi, Mi m and A;o must all be considered as independent. The corresponding constrained (representation of 
the) instant-form generators, with po and A;o expressed in terms of the constrained state x', become therefore 



Pol 



Nu 



E 

i=l.N 

E 



/ (»)2 2 

' m n c z 



c (e//) 



g W 4(«ot)(i) 
c (e//)0 



=1JV 



(») / (02 2 

r i A / m ° 



PIM _ A (fot)(iA , 9 W (i) A(tot)(i) 



(153) 
(154) 
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where we have represented A*f°fl\ 1 ^ = C^^ff)o> A*'* ' ''" 



with A 



{totm being defined by Eq.ftD setting A^ xt)[i) = 0. 



~ \ (eff)O^(eff) J' ^(eff)n 

A characteristic obvious feature of the constrained representations given above is that of the non-local dependences 
arising both from binary and self EM interactions. Analogous conclusions can be drawn also for the so-called point 
and front forms of the same generators. This implies that the Lorentz conditions (|126|) . with PBs now defined in 



in this case 



(155) 



terms of the same constrained state x', are generally violated. Indeed, due to the non- locality of ^[^//j/* 
the PBs-inequalities 

[pox',p;] (x0 + 

hold. Hence, if - consistent with DGF - the validity of the Lorentz conditions \126\) is imposed, the constrained forms 
of the Poincare generators are manifestly not applicable to the treatment of EM-interacting N-body systems. 

Nevertheless, it is immediate to prove that po\ x , indeed generates the correct evolution equations for the constrained 
state x'. In fact, denoting by [•, -b x ,) the local PBs evaluated with respect to the constrained state x', let us determine 
by means of the PBs 



«U(x) = [£,F] 



(x') 



the infinitesimal transformations 5 r^ and 5 P^ generated by F = dtpo\ 
yield respectively 



(156) 

It is immediate to prove that these 



<5 r«^ = dt 
<5 D P« = dt P«,p 



(X') 



EE (ftyW, 



(x') 



(157) 
(158) 



where the r.h.s. of both equations coincide identically with the spatial parts of the canonical Eqs. (j9"?| and ([95]) . Hence, 
as expected, the constrained state x' is indeed canonical. Eqs. |-?57[ ) and \158\) provide the Hamiltonian equations for 
x' in terms of the non-local Hamiltonian function pb| x /- Again, a necessary and sufficient condition for Eqs. (|97p and 
to hold is that the fundamental PBs 



r («) ; r (i) 




= o, 




K) 




p(») p(j) 




= o, 




K) 


r (*) pO') 


K) 


= S ij l, 



(159) 

are identically satisfied for all i,j = 1,N. Here Xq and x' are identified respectively with x' = x'(t Q ) and x' = 
x'(t) = x' a + dx', with dx'= (i^oX^ 1 ', .., S a yS N ^, while the previous PBs are evaluated with respect to the initial state 
xj,. Furthermore, also in this case the canonical initial history set {x'} to , to be defined in analogy with Eq. ([M|) . is 
assumed prescribed. 



IX. NON-LOCAL GENERATOR FORMALISM 



A basic consequence of the previous considerations is that in the case of non-local phase-functions, such as Hn or 
Po\ x r, the local transformation law (| 1 29[1 becomes inapplicable. 

A suitably-modified formulation of DGF appropriate for the treatment of non-local phase-functions must therefore 
be developed. This can be immediately obtained. In fact, let us consider an arbitrary non-local function of the form 
£ = £(z, [z]), with z and [z] denoting respectively local and non-local functional dependences with respect to the 
canonical state z. Let us consider an arbitrary infinitesimal canonical transformation generated by F of the form 
z — > z' = z + [z, F], with 5 z to be considered as infinitesimal (i.e., of 0(A)). Then, requiring that £ is suitably 
smooth both with respect to z and [z] , the corresponding infinitesimal variation of £ can be approximated with 

<^(z, [z]) ee [£(z + a6 z, [z + a5 z]) - £(z, [z])] [1 + 0(A)] , (160) 

6£(z, [z]) being the (Frechet) functional derivative of £(z, [z]), namely 

<S£(z, [z]) ee lim (z + aS z, [z + aS z]) = {£(z, [z]), F} . (161) 
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Here {£(z, [z]), F} denotes the non-local Poisson brackets (NL-PBs) and generally also F can be considered a non-local 
function of the form F (z, [z]) [i.e., of a type analogous to £]. Such a definition reduces manifestly to (|129[) in case of 
local functions. 

Let us prove that the transformation law (I161[) is indeed the correct one. To elucidate this point, let us consider the 
4-scalar defined by the Dirac-delta £(r, [r]) = S (lt^ a Ra^ — °fi)J entering the non-local Lagrangian and Hamiltonian 

functions in the EM self- interaction, where R^ a denotes the bi- vector defined by Eq. (f2"5)l . Let us consider, for 

example, the action of an arbitrary infinitesimal Lorentz transformation defined by 8 r^^ . In order that PS l ^ a PSa is 
left invariant by the transformation (Lorentz invariance) it must be 

tf£(r,[r]) = {£(r ) M),F} = ) (162) 

with F = -p M a M . This means that the NL-PBs {£(r, [r]), F} defined by Eq. pMj) . rather than the local PBs 
[£(r, [r] ) , F] , must vanish identically. In particular Eq. (116ip , contrary to the local variation (|128[) , preserves the 
Lorentz invariance of 4-scalars and hence provides the correct transformations law. Hence, in particular, it follows 
that for an isolated TV-body system with arbitrary N > 1 : 

8H N (r,P,[r}) = {H N (r, P, [r]), F} — 0, (163) 
dp \ x , = {po\ x ,,F} = 0. (164) 

It is immediate to prove that Eq. (|163[) holds by construction for all Poincare generators [see Eqs. (|148p above], while 
- instead - generally 

S H N (r,P,[r}) ee [H N (r,P,[r]),F}^0, (165) 
5oPo\ x , = [PoU^^O. (166) 

Hence, consistent with the results indicated above, we conclude that the local transformation laws realized by the 
Lorentz conditions U 2b}) . which are a distinctive feature of DGF, become invalid in the case of non-local Hamiltonians. 

The non-local generator formalism is therefore formally achieved by imposing modified Lorentz conditions obtained 
from Eqs. (|126l) . in which the local PBs are replaced with the non-local PBs defined by Eg. M6l\ ). 

In particular, the correct transformation laws for the constrained instant-form Poincare generators [see 
Eqs. (|150l) . (11511) and (|153[) . (|154[) ] follow by imposing for the Hamiltonians [see Eqs. (ll53[) and (|154p ] appropriate non- 
local transformation laws of the type \16J$ , all defined with respect to the constrained state x'. Finally, it must be 
remarked that the non-local generator formalism does not affect the validity of the canonical equations of motion 
(|157|) and (|158|) as well as the fundamental PBs (|159|) indicated above for the constrained state x', which remain 
unchanged. 



X. COUNTER-EXAMPLES TO THE "NO-INTERACTION" THEOREM 

An open problem in relativistic dynamics is related to the so-called "no-interaction" theorem due to Currie (Currie, 
1963 jJL3(), derived by adopting the DGF, and in particular the instant form representation for the Poincare generators 
(see previous Section) given in Ref.[3]. According to this theorem, an isolated classical A-body system of mutually 
interacting particles which admits a Hamiltonian structure in which the coordinate variables of the individual particles 
coincide with the space parts 3-vectors of the particles 4-positions and the canonical equations of motion are Lorentz 
covariant, can only be realized by means of a collection of free particles. This requires, in particular that "...it is 
impossible to set up a canonical theory of two interacting particles in which the individual particle positions are the 
space parts of ^-vectors" . In other words, according to the theorem, it should be impossible to formulate - in terms 
of a Hamiltonian system - a covariant canonical theory for an isolated system of N > 1 classical particles subject to 
binary interactions (see also Ref . fl8l ) . The validity of the theorem was confirmed by several other authors (see for 
exam ple, Beard and Fong, 1969 [23[, Kracklauer, 1976 (2^, Martin and Sanz, 1978 [25], Mukunda and Sudarshan, 
1981 26], Balachandran et al, 1982 [27]). Its original formulation obtained by Currie for the case of two interacting 
particles (N = 2) was subsequently extended to include the case A = 3 (Cannon and Jordan, 1964 [28|), first-class 
constraints (see Sudarshan and Mukunda, 1983 [2!| and the corresponding Lagrangian proof given by Marmo et al., 
1984 (30[) and the treatment of curved space-time (De Bievre, 1986 [3lJ and Li, 1989 (32j). Common assumptions to 
these approaches are that: 
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1. Hypothesis #1: Both DGF and the Dirac instant form realization of the Poincare generators apply. In particular, 
the Poincare generators in the instant form, corresponding to the constrained Hamiltonian structure {z, JiTjv} , 
satisfy identically both to the commutation rules (|126|) and the kinematic constraints (|10ip . 

2. Hypothesis #2: Kn admits the Poincare group of symmetry, i.e., it commutes with {F}. 

3. Hypothesis #3: All particles, in a suitable proper-time interval, are not subject to the action of an external 
force (locally/globally isolated TV-body system). 

Nevertheless, the theorem has been long questioned (see for example Fronsdal, 1971 |l5| and Komar, 1978-1979 
[l6l - [l9l ] ) . In particular, there remains the dilemma whether the "no- interaction" theorem actually applies at all for 
TV-body systems subject only to non-local EM interactions. This refers in particular, to extended charged particles 
in the presence of binary and self EM forces. Another interesting question is whether restrictions placed by the 
"no-interaction" theorem actually exist for physically realizable classical systems. Several authors have advanced the 
conjecture that the limitations set by the Currie theorem might be avoided in the framework of constrained dynamics 
formulated adopting a super-abundant-variable canonical approach (see for example Komar, 1978 [l8T ] and Marmo et 
at, 1984 [3fJ and references indicated therein). In particular, to get a better understanding of interacting A-body 
systems, Todorov [33| and then Komar |16l - [l9l | developed a manifestly covariant classical relativistic model for two 
particles, of an action-at-a-distance kind. In the Todorov-Komar model the dynamics is given in terms of two first- 
class constraints. An equivalent model was discovered by Droz- Vincent [34j,[35j based on a two-time formulation of the 
classical relativistic dynamics. However, the precise identification of the Hamiltonian structure {z,K^} pertaining to 
A-body systems subject to EM interactions has remained elusive to date. 

Here we claim that counter-examples, escaping both the assumptions and the restrictions of the "no-interaction" 
theorem, can be achieved, based on the classical A-body system of extended charged particles formulated here. 
Starting from the Corollary to THM.2, the following theorem applies. 

THM.4 - Standard Hamiltonian form of a locally-isolated 1-body system and a globally-isolated 
A-body system. 

In validity of THM.2 and of the definitions given by Eqs. |7_?| )- J76p , the following propositions hold: 

T4i) The Hamiltonian structure {x, iJ/v } of the classical system formed by a single extended charged particle is 

preserved also in the particular case in which the external EM A—potential is such that along the particle world-line 

r«( S(1) ); 

4 ex4) (r- (1) ( S(1 ))) = (^ n ° Vs M e ]-°o, So ] (167) 

M K > \ vs(i) e \s , +00, [ v ' 

(locally-isolated particle). 

TJj.<z) The Hamiltonian structure {x, if/v} of the classical N-body system formed by extended charged particles is 
preserved also in the particular case in which the external EM A—potential vanishes identically, 

A^ xt \r) = (168) 

(globally-isolated N-body system). 
Proof - T4i) The proof is an immediate consequence of the Corollary to THM.2. In fact in the absence of an 

external EM field, the effective EM 4-potential A?°fj^ (see Ea.([73]l) simply reduces to 

4°mJ = 2 t eim \ (169) 

where, in view of the requirement (|167|) . ^ ^ is non- vanishing also in the interval ]s , +00, [ (see related discussion 
in Paper I). Hence, both the Lagrangian and Hamiltonian equations in standard form [see respectively Eqs. (|57[) and 

flM]),®] are satisfied, with H^ f = #jv,e// = Hn still defined by Eqs.flTf)]) and (JMJ). T4 2 ) The proof is similar. In 

this case, due to assumption f|168[) . A^°Jj^ reduces to 

(i) 

Hence, also in this case both the Lagrangian and Hamiltonian equations in standard form still hold, with H e ^ and 
HN,eff defined by Eqs. (US]) and 
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Q.E.D. 

It is clear that both propositions T4i) and T4 2 ) indeed escape the "no interaction" theorem (avoiding also the 
limitations set by its assumptions #l-#4). In fact, concerning the Hamiltonian structure {x, Hn} associated to the 
classical iV-body system of extended charged particles, from THM.4 it follows that: 

• The effective Hamiltonian is a non-local function of the canonical state x. 

• The canonical particle equations of motion (I88[) satisfy the correct transformation laws with respect to the 
Poincare group, since the non-local system Hamiltonian Hjy(r, P, [r]) is by construction a Lorentz 4-scalar. 

• The super-abundant canonical state x = {x^,i = 1, AT} is defined in terms of = (r^^^P^j , where 

v ' ( s (i)) 

rW'* and Pji' are represented by Eqs. (|90| and (|9l"j) . 

• The extremant curves (x' 1 ' (sm), ...x^ N ' (s(jv))) solutions of Eqs. (|123|) satisfy identically the kinematic con- 
straints (jlOip and (|102p . As a consequence, only the space parts of the extremant 4-vectors r"''(sn)) and 
Pff 1 (s(i)) are, for all i = l,N, actually independent. 

• In case of T4i), the single-particle motion is non-inertial for all sm £ / = R. Hence, the instant form of Dirac 
generators [see Eqs. (ll34[ ) and (|137[) ] becomes inapplicable even in the case of a 1-body system. 



XI. ON THE FAILURE OF THE "NO-INTERACTION" THEOREM 



The actual causes of the failure of the "no-interaction" theorem emerge clearly from the analysis of the conditions 
of validity of DGF and the Dirac instant-form generators (see Section 8). For systems of extended charged particles 
subject only to EM interactions the previous assumptions #l-#4 (common to all customary approaches [TH, [23l — 
l3l| ) which characterize the underlying Hamiltonian structure {z, Kjy} make it incompatible with the exact non-local 
Hamiltonian structure {x, Hn} determined here. In fact, in difference to {z, Kn}, the Hamiltonian structure {x, H^} 
is characterized by: 

• Super-abundant canonical variables x = \x.^\i = l,N\ , with x^ = (r^^P^lr ) being the i-th particle 
canonical state. 

• Extremant curves (xW^m), ■•• x ^( s w)) which satisfy identically the kinematic constraints discussed above. 
This is a characteristic property of the canonical extremant curves only. In fact, the same constraints are not 
satisfied by the super-abundant canonical state x. 

• Fundamental PBs (|104l) which are satisfied only by the unconstrained state x. Hence, the non-local Hamiltonian 
structure {x, Hn } is warranted if all the canonical variables defining the state x = {x^- 1 ,i = 1, A'} are considered 
independent. This means that, in order for the fundamental PBs (|104l) to be fulfilled, these constraints cannot 
be imposed "a priori" on the canonical state. 

• Poincare generators [see Eqs. (ll49j) - (ll52p ] which satisfy the commutation rules (|126l) when they are considered 
independent, as the super-abundant canonical variables x = {xW,i = 1, AT} , and fulfilling the fundamental 
PBs. For this reason, the Poincare generators are necessarily left unconstrained by imposing the validity of the 
same equations [i.e., Eqs. ([126[) ]. 

• A non-local Hamiltonian of the form -ff/v(r, P, [r]). In particular, it follows that Hn for classical A^-body systems 
of extended charged particles reduces to a local function only in the case of a single isolated particle which exhibits 
inertial motion. In view of THM.l given in Paper I, this requires the external EM 4-potential acting on such a 
particle to vanish identically along the particle world-line, i.e., A^ xt \r{si)) = for all s\ £ / = R. 

In conclusion, contrary to the claim of the "no-interaction" theorem, a Lorentz covariant Hamiltonian formulation 
for the dynamics of N-body systems, with N > 1, actually exists also for mutually interacting charged particles subject 
to binary as well as self EM interactions. The result holds even in the presence of an external EM field, for extended 
classical particles described by the Hamiltonian structure {x, Hn } determined here. 
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One might conjecture that the validity of the "no-interaction" theorem could be restored by introducing a suitable 
asymptotic approximation for the A*-body system dynamics. The latter is related, in particular, to the short delay- 
time and large-distance approximations (see Section 7) , invoked here for the treatment of particle self and binary EM 
interactions. It is immediate to prove that also this route is necessarily unsuccessful. The reason lays in THM.3 and 
its consequences. In fact, as shown above, a Hamiltonian structure of the same type of {x, Hn} can be recovered 
for the asymptotic A^-body equations of motion determined by the same theorem. This is identified with the set 

{*> H NTf)}> where H N V e?f = Ei=i,jv H ef /.cwym and H eff,as ym is § iven h Y Eq.<H2D>- By construction, {x,H£%fj } 
inherits the same qualitative properties of the exact Hamiltonian structure {x, H^}. Therefore, in particular, in this 
approximation x satisfies the fundamental PBs (|104j) if it is unconstrained. In addition, since the same definition 
applies for the Poincare generators and their representation in the instant form, the same conclusions on the validity 
of the "no-interaction" theorem follow. 



XII. CONCLUSIONS 

A formidable open problem in classical mechanics is provided by the missing consistent Hamiltonian formulation 
for the dynamics of EM-interacting Af-body systems. This critically affects both classical and quantum mechanics. 
In this paper a solution to this fundamental issue has been reached exclusively within the framework of classical 
electrodynamics and special relativity. In particular, the Hamiltonian structure of classical A^-body systems composed 
of EM-interacting finite-size charged particles has been explicitly determined and investigated. 

Both local and non-local EM interactions have been retained. The former are due to externally-prescribed EM fields, 
while the latter include both binary and self EM interactions, both characterized by finite delay-time effects. Binary 
interactions occur between any two charges of the A"-body systems, while self interactions ascribe to the so-called 
radiation-reaction phenomena due to action of the EM self-field on a finite-size particle. All of these contributions 
have been consistently dealt with in the derivation of the Af-body dynamical equations of motion by means of a 
variational approach based on the hybrid synchronous Hamilton variational principle. 

Both Lagrangian and Hamiltonian covariant differential equations have been obtained, which are intrinsically of 
delay-type. The same equations have also been proved to admit a representation in both standard Lagrangian and 
Hamiltonian forms, through the definition of effective non-local Lagrangian and Hamiltonian functions. The property 
of Hamilton equations of admitting a Poisson bracket representation has lead us to prove the existence of a non- 
local Hamiltonian structure {x, Hjy} for the A^-body system of EM-interacting particles. This has been shown to 
be determined by the non-local Hamiltonian function iJ/v and to hold for the superabundant canonical states x. In 
particular the correct Hamiltonian equations of motion are obtained considering the same vector x as unconstrained, 
the relevant (kinematic) constraints being satisfied identically by the solution of the same equations. 

A further interesting development concerns the asymptotic approximation determined for the Hamiltonian structure 
{x, iJjv} of the full A^-body problem. Here we have shown that consistent with the short delay-time and large- 
distance asymptotic orderings the latter can be preserved also by a suitable asymptotic Hamiltonian approximation. 
In particular, the perturbative expansion adopted here permits to retain consistently delay-time contributions, while 
preserving also the variational character and the standard Lagrangian and Hamiltonian forms of the A^-body dynamical 
equations. As a basic consequence the very Hamiltonian structure of the A^-body problem is warranted. This permits 
us to overcome the usual difficulties related to the adoption of non- variational and non-Hamiltonian approximations 
previously developed in the literature. 

Two important applications of the theory have been pointed out. 

The first one concerns the famous and widely cited (both in the context of classical and quantum mechanics) paper 
by Dirac (1949) on the generator formalism approach to the forms of the Poincare generators for the inhomogeneous 
Lorentz group. Contrary to a widespread belief, we have found out that the Dirac approach is not valid in the case 
of A^-body systems subject to retarded, i.e., non-local, interactions. In fact, the Lorentz conditions for the instant- 
form Poincare generators are found to be satisfied only in the case of local Hamiltonians. Analogous conclusions can 
be drawn also for the so-called point and front-forms of the same generators. Due to the non-local character of the 
Hamiltonian structure {x, H^} this means that the Dirac generator formalism expressed in terms of the essential (i.e., 
constrained) canonical state x' is not "per se" directly applicable to the treatment of EM-interacting A^-body systems. 
However, as shown here, in the same variables its extension to non-local Hamiltonians can be readily achieved by 
suitably modifying the Lorentz conditions so to account for the non-local dependences of the Hamiltonian structure 
{x,H N }. 

Second, the validity of the Currie "no-interaction" theorem, concerning the Hamiltonian description of the relativis- 
tic dynamics of isolated interacting particles, has been investigated. It has been proved that the set {x, Hn} violates 
the statements of the theorem. The cause of the failure of theorem (and its proof) lays precisely in the adoption of the 
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Dirac generator formalism. Explicit counter-examples which overcome the limitations posed by the "no-interaction" 
theorem have been issued. Contrary to the claim of the "no-interaction" theorem, it has been demonstrated that 
a standard Hamiltonian formulation for the iV-body system of charged particles subject to EM interactions can be 
consistently formulated. 
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Appendix A: Evaluation of the action integral of the binary interaction 



In this Appendix the mathematical details of the calculation of the action integral 5p m ^^(r, [r]) are given. The 
latter contains the information about the binary EM interactions among the charged particles of the iV-body system 
and it has been defined by Eqs. (fT4"| and (fl"5|) in Section 4- A. To proceed with the calculation we first notice that, 
invoking the definition of the current density given by Eq.(j3]), the functional (|15p can be equivalently represented as 



S, 



(bin)(ij) 
C 



(r, [r]) 



dVtA (seim ^ l (r) 



ds 2 S(s 2 



5 i(i)J 



+ 00 



dsyjuW" (s U) )6(\x (s u) ) I - o- U) )6{s (j) 



S2), 



where Suj) is the root of the equation 



uj^faiCj)) 



= o. 



(Al) 



(A2) 



Furthermore, because of the principle of relativity, the integral (|Alj) can be evaluated in an arbitrary reference 
frame. The explicit calculation of the integral (|A1[) is then achieved, thanks to Lemma 3 given in Paper I, by 
invoking a Lorentz boost to the reference frame Sni moving with 4- velocity u^(s 2 ). In this frame, by construction 
dCl' = cdt' dx' dy 1 dz' = dft. In particular, introducing the spherical spatial coordinates (ct' , p', <p') it follows that the 
transformed spatial volume element can also be written as cdt' dx' dy' dz' = cdt' dp' dd' dip' p' 2 sin??'. In such a reference 
frame the previous scalar equation becomes 



= 0. 



(A3) 



On the other hand, performing the integration with respect to s 2 in Eq. (|Aip , it follows that necessarily s 2 — Si(j), so 
that from Eq. (|A3p Si(j) is actually given by 



ct' 



s 2 - 



As a result, the integral Sq"' ^ reduces to 



S ( c nKij) (r',[r'}) 



-,(./) 



47TCT 



(i) 1 



+oo 



dx'dy'dz'A( sel fW»( r >) j ds {j) u^ ,fl (s {j) )S( x&' ( S(j) ) 

with xW'i* (s 0) ) = r'^ - r^ 1 (s {j) ). Moreover 

A ^m) V) = 2g W / ' ™ ds '^V (s'fo) 5 {R^R^') 



(A4) 



(A5) 



(A6) 
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with R^' a = r' a - r^' a (s'',) 



(i)> 

Hence, S^ m ^ 1 ^ reduces to the functional form: 



q q 1 dtf'sintf' 



or 



S^WJ^.tr 7 ]) = / dtf'sintf' / V / dp'p' 2 x 



+ 00 



ds^f (V ( '. } ) [ + °°ds {j) u^( S{j) )5(xW{s {j) ) -a {j) ). (A7) 

J —00 



The remaining spatial integration can now be performed letting 

p'=\x' (s (j) )\ 



(A8) 



and making use of the spherical symmetry of the charge distribution. The constraints placed by the two Dirac-delta 
functions 5{R^' a R^a') and 5{\x' (s(j)) | — in the previous equation imply that both R^' a Ra^' and \x' | are 
4-scalars. Then, introducing the representation 



£(i)/a = r /a _ r W'«( s ' ( '. ) ) = fl(«)'a + ^/a ^ ^ 



(A9) 



with 



r(W>* = r a)'« ( Sa) ) _ r W'«( s ' ( '. ) ) ; 



(A10) 

= r~-r»'~(8 U) ), (All) 
it follows that 

R(i)'afi(i)t = R«i)'*R&)> + X U)»* ( S(j) ) X U)> ( S(j) ) + 2R^ /a x^' (s {j) ) (A12) 

is necessarily a 4-scalar independent of the integration angles (<^',i?') when evaluated on the hypersurface E : 

= _ Similarly, the Dirac-delta (s (j) )| - cr 0) ) warrants that (s (j) ) (s (j) ) 

is manifestly a 4-scalar too. Let us now prove that necessarily 

rW)'<* x U)> ( S(j) ) = 0. 

In fact, on E it must be 



-a 2 {j y which 



d 



dsa 



R(i)"*R(i)> 



ds" 



0. 



ds 



ds" 

as (i) 



i^>^>'( S(j) ) 
d 



u^' a (s {j) )x^'(s {j) )-R^u^'(s {j) ) 
1 d 



2ds 



0) 



R(ij)la R(ij)< 



ds" 



R(ij)'acR(ij)l 



Therefore, 



d 




d r 






ds (j) L 



-« (i),a (4)) ^' (-W)) . 



rW)"*rW)' + 2R (ll),a x^' (s U) ) 



(A13) 
(A14) 

(A15) 
(A16) 
(A17) 

(A18) 



ds" 



R(i)ICCR(i)> 



ds" 

as (i) 



rW)'<*rW)' + 2R (l ^' a x ( ^' (s 0) ) 



-2fl«> to ««' (V ( ' i} ) - 2 U W- (^) 



= 0, 



(A19) 
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from which it follows that, on £, is a 4- vector, since by definition both u^' (s"^ and (s (j) ) are 4- vectors 

too. Now we notice that 

Rm»>RW = / ( S(J) , S ' ( ' ) = / ( S ' ( ' i} , S {j) ) , (A20) 

with / being a 4-scalar which is symmetric with respect to sy) and s"^ , while by construction 

RM"*x$'{s m ) =g[s {j) ,s'l {) ,a {j) ) ^ g (s'fo,s {j) ,a {i) ) , (A21) 

where g is a non-symmetric 4-scalar with respect to the same parameters. On the other hand, Eq. (|A14[) requires that 
R^' a Ra^ must be symmetric in both S(j) and s'/^, so that, thanks to Eqs. (|A20p and (|A21[) . we can conclude that g 
is a constant 4-scalar. To determine the precise value of g — R^' a Xa ( s (j)) we evaluate it in the j-th particle COS 
co-moving reference frame, where by definition tqq S (sq) = (so, 0) for all the COS proper times sq 6 [— oo, +oo]. In 
this framei?^*"" = \S(j^ — s"^,0^ has only time component and when sq — we get g = R^' a x^' (sy)) = 
identically. On the other hand, since g is a 4-scalar, it is independent of both sy) and s"^ and it is null when so = sy), 
we conclude that it must be null for all so and in any reference frame, which proves Eq. (|A13j) . 

Hence, as a result of the integration, the action integral Sq 71 ^ 1 ^ carrying the interaction of particle i on particle j 
takes necessarily the expression 

S^ ij \r', [/]) = fdrf f^r^is^m^R^' - (A22) 



Finally, since by construction S^" 1 ^ 1 ^ is a 4-scalar, by dropping the primes and replacing s'/* and s',^ respectively 
with S(j) and srj\, the result reported in Eq. (|16[) is recovered. 
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